Spin transfer in bilayer magnetic nanopillars at high fields 
as a function of free-layer thickness 
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Spin transfer in asymmetric Co/Cu/Co bilayer magnetic nanopillars junctions has been studied 
at low temperature as a function of free-layer thickness. The phase diagram for current-induced 
magnetic excitations has been determined for magnetic fields up to 7.5 T applied perpendicular to 
the junction surface and free-layers thicknesses from 2 to 5 nm. The junction magnetoresistance 
is independent of thickness. The critical current for magnetic excitations decreases linearly with 
decreasing free-layer thickness, but extrapolates to a finite critical current in the limit of zero 
thickness. The limiting current is in quantitative agreement with that expected due to a spin- 
pumping contribution to the magnetization damping. It may also be indicative of a decrease in the 
spin-transfer torque efficiency in ultrathin magnetic layers. 

PACS numbers: 



Spin transfer in magnetic nanopillar has become a ma- 
jor focus of experimental research [J 0, 9 since Slon- 
czewski and Berger's seminal theoretical work in 1996 
[1,0. A spin current has been demonstrated to switch 
the magnetization direction of a small magnet at a spe- 
cific current density, as well as to induce microwave exci- 
tations. There are applications of this effect to magnetic 
random access memory (MRAM) and high-frequency 
electronics [H @> 0- ^ is of importance to determine 
the factors that control the critical current for magneti- 
zation dynamics for both the physics and technology of 
spin transfer. For instance, it is of interest to reduce the 
critical current for MRAM applications, and to increase 
it in magnetic sensor designs. 

In Slonczewski's theory, spin transfer is an interface ef- 
fect: spin- angular momentum is transferred to the back- 
ground magnetization when the spin current enters the 
ferromagnet -within the first few atomic layers Q. For 
one polarity of the current, this generates a torque on the 
magnetization that is opposed by bulk damping. As a re- 
sult, there is a threshold current to excite magnetization 
dynamics that is proportional to the volume of the mag- 
net or, equivalently, the threshold current density is pro- 
portional to the thickness of the magnetic layer. There 
are alternative models in which the spin-transfer interac- 
tion occurs on a longer length scale, which predict a de- 
crease in the efficiency of the torque in very thin magnetic 
layers [j| [l(| • It is also now widely appreciated that the 
magnetization damping in thin layers can be dominated 
by interfaces, in an effect known as "spin pumping" [Til ] . 
For these reasons it is of importance to study spin trans- 
fer in samples in which the layer thicknesses are varied to 
gain insight into the factors that determine the strength 
and length scales of the spin-transfer interaction. 

Albert et al. _12j studied current-induced magnetiza- 
tion switching as a function of free-layer thickness at 
room temperature. Here thermal fluctuations are im- 
portant and the intrinsic (zero temperature) critical cur- 
rent was determined by extrapolating from pulsed cur- 



rent measurements. The switching was between in-plane 
magnetized states, parallel and antiparallel to the fixed- 
layer magnetization. In this case, the in-plane shape 
anisotropy plays an important role in setting the energy 
barrier to reversal. The switching current was found to 
depend linearly on the free-layer thickness and to be zero 
in the limit of zero free-layer thickness, consistent with 
Slonczewski's model. 

In this paper, we present studies of spin-transfer at 
low temperature and high magnetic fields in asymmet- 
ric magnetic nanopillars in which the free magnetic-layer 
thickness has been systematically varied. The phase di- 
agram for magnetic excitations has been determined in 
fields perpendicular to the film plane, under which the 
in-plane anisotropy is a minor effect, and at low tem- 
perature (4.2 K), where thermal fluctuations can be ne- 
glected. Similar to the results of Albert et al. [Hj], we 
find that the critical current density, defined to be the 
current density at which there is a step change in junc- 
tion resistance, is a linear function of the layer thickness. 
However, it extrapolates to a finite zero thickness inter- 
cept. This suggests that damping related to spin pump- 
ing sets a lower limit for the critical current density in 
ultra-thin magnetic layers, or that spin transfer occurs 
over a finite-length scale in the ferromagnet. 

In magnetic nanopillars that consist of thick and thin 
magnetic layers separated by a nonmagnetic layer, the 
thick (fixed) layer polarizes the current and dynamics is 
induced in the thin (free) layer. The magnetization of 
the free layer can be described by the Landau-Lifshitz- 
Gilbert (LLG) equation with an additional spin-transfer 
torque term [J, |5(. In a macrospin model, which assumes 
that the two layers are uniformly magnetized: 
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m and m p are unit vectors in the direction of magnetiza- 
tion of the free and fixed magnetic layers, respectively. 7 
is the gyromagnetic ratio. The second term on the right 
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Jt~ap = ^p-(l - v)M s t(H + 4vrM cff ). 
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FIG. 1: Phase diagram of a nanopillar for 77 = 0.3, a = 0.01, 
and 47rM c ft = 1.5 T, with a small (0.1 T) in-plane uniaxial 
anisotropy. The threshold current densities for instability of 
the P state, Jt-p, and for instability of the AP state, Jt-ap, 
are indicated. 



is the damping term, where a is the Gilbert damping 
constant. The last term is due to spin-transfer. Here 
a j = 2eM°t ' wnere J is the current density, P is the 



spin-polarization of the current, M s is the magnetiza- 
tion density, and t is the free magnetic-layer thickness. 
The positive constant r\ characterizes the angular depen- 
dence of the torque and depends on the spin polariza- 
tion [lj| . The effective magnetic field H e g is the vec- 
tor sum of the applied held H, the in-plane anisotropy 
field -f/a, and the easy-plane anisotropy —&TrM e ffm z z, 
H c g = H + H a m x x — AnM c tfm z z. In equilibrium, the 
magnetization is aligned with the effective magnetic field. 
The spin-transfer torque competes with damping, and 
at a threshold current density leads to excitation of the 
magnetization. 

A phase diagram assuming a single-domain macrospin 
model can be calculated from Eq. 1, as shown in Fig. [1] 
For large fields perpendicular to the plane of the magnetic 
layers (H > AnM s ), the magnetizations of the fixed and 
free layer are aligned in the field direction. This is a 
particularly simple situation [2J, as the in-plane shape 
anisotropy of the element plays only a minor role in the 
dynamics -the easy-plane anisotropy dominates, and the 
magnetic energy has axial symmetry. In this case, as J 
increases, the parallel (P) state becomes unstable at a 
threshold current: 
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leading to a precessional (PS) state. A further increase 
in the current results in the free layer switching into an 
antiparallel (AP) state. For a decreasing current, the AP 
becomes unstable and switches back to a PS or P state 



Jt-p and Jt-ap cross and become equal at the field 
inM e ff/r]. However, hysteresis appears at a smaller field, 
H* , corresponding to the lowest field at which the AP/PS 
region appears in the phase diagram. An interesting fea- 
ture of the perpendicular field phase diagram is that the 
hysteresis is associated with the angular dependence of 
the spin-transfer torque the fact that the torque is 
larger in the AP state (77 > 0). This is in contrast to 
conventional hysteresis in magnets, which is associated 
with dipolar interactions or magnetic anisotropy. Note 
that starting at large current and decreasing the current, 
the model predicts an abrupt and large change in the an- 
gle between the free and fixed layer at Jt-ap, which is 
detectable as a step change in junction resistance. This 
suggests that a measurement of the critical current on 
a decreasing current can be used to determine the AP 
threshold current (Eq. 3) as a function of magnetic field 
and sample structure. This is the approach we take in 
analyzing our experiments. While we will focus on the 
high-field behavior, we note for fields less than 4ttM c s the 
layer magnetizations tilt into the film plane and the phase 
boundaries depend on the in-plane magnetic anisotropy. 

Hundreds of pillar junctions with submicron lateral di- 
mension were fabricated on a 1 x 1 cm 2 silicon wafer using 
a nanostencil process [lij]. Stencil holes with different 
but accurate lateral dimensions were opened up at the 
depth of ~ 75 ran, and pillar junctions were deposited 
into those stencil holes through metal evaporation. Junc- 
tions have the layer structure || 3 nm Pt | 10 nm Cu \ 
t Co I 10 nm Cu | 12 nm Co j 10 nm Cu | 3 nm Pt j 
200 nm Cu ||. During the evaporation of the thin Co 
layer, we used a linear motion shutter to vary t from 1.8 
to 5.3 nm across the wafer. Junctions with lateral dimen- 
sions 50 x 50, 50 x 100 and 70 x 140 nm 2 were studied 
in detail. 

All measurements reported here were made at 4.2 K 
with a 4-point geometry. Both dc resistance V/I and dif- 
ferential resistance dV/dl were measured for each junc- 
tion. A 0.2 mA modulation current at 802 Hz was added 
to the dc bias. Junction resistances were found to scale 
inversely with lateral areas. Positive current is defined 
to be electron flow from the free (thin) Co layer to the 
fixed (thick) Co layer. 

The magnetoresistance (MR) was measured with the 
magnetic field applied in the film plane. A typical MR 
hysteresis loop of a 50 x 50 nm 2 junction with t ~ 2.8 nm 
is shown in Fig. [2ja) . The high resistance state corre- 
sponds to an AP state and for fields greater than 0.1 T 
lead to a P state of lower resistance. The magnetoresis- 
tance, MR = {R AP -Rp)/R P = SR/Rp, is (2.2±0.2)%, 
independent of the free layer thickness within the thick- 
ness range investigated (Fig. 1(c)). This shows that the 
MR is due mainly to spin-dependent scattering at Co/Cu 
interfaces. The MR area product, SRA, where A is the 
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FIG. 2: (a) Zero dc current in-plane magnetoresistance hys- 
teresis loop for a 50 x 50 nm 2 junction with t ~ 2.8 nm. (b) 
Positive current sweep hysteresis loop of the same junction 
with perpendicular magnetic field set at 7 T. (c) MR of all 
junctions as the function of t. (d) Zero dc current in-plane 
SR times lateral area A for all junctions as the function of t. 



lateral area of the junction, is also independent of thick- 
ness, which indicates that SR is inversely proportional to 
the junction expected. 

Current-voltage measurements were conducted with a 
magnetic field applied perpendicular to the sample sur- 
face. Fig. [2^b) shows a current sweep hysteresis loop of 
the same junction as in Fig. [2a) with a 7 T applied field. 
When the current is swept up to a sufficiently large value, 
17.1 mA in this case, a step increase in resistance is ob- 
served (indicated with an upward arrow in Fig. [2b)). 
We interpret this step as indicating the current at which 
the junction switches into an AP state, as reported in 
Ref. 0. For decreasing current there is a step down in 
differential resistance at 15.7 mA, which, as discussed, 
we associate with the linear instability threshold given in 
Eq. [31 Jt-ap- The current hysteresis is about 1.4 mA 
at this field. The majority of junctions show hysteresis 
in current sweep measurements for fields larger than 3 T 
and steps in both dV/dl and V/I 15]. The sharp step in 
V/I has been used to determine the critical current for 
all junctions. 

A contour plot of dV/ dl as a function of current and 
perpendicular magnetic field shows the variation of the 
critical current with the applied perpendicular field. In 
order to emphasize the change of resistance on top of 
the background, which is associated with Joule heating, 
we have plotted dV/dl minus a linear background ver- 
sus current density. Fig. [3£a) shows data for decreasing 
current on the same junction as in Fig (2a), (b). The 
brighter the color in the contour plot, the larger the junc- 
tion differential resistance. The boundary between the 
bright and dark region is an abrupt step in differential 
resistance. The corresponding position of the abrupt step 
for increasing current is illustrated with black dots in Fig. 



FIG. 3: (a) Contour plot of dV/dl minus a linear background 
as a function of both current density and magnetic field per- 
pendicular to sample surface for decreasing current. The same 
junction as in Fig. 2 (a,b). "•" data points: The correspond- 
ing step in V/I for increasing current, (b) V/I vs current bias 
hysteresis loops with fields set at 3, 5 and 7 T. (c) SR/R at 
the critical current as the function of field. Solid line: Zero 
dc current in-plane MR. 



EJa). We note that this boundary is not as straight as 
the current sweep down boundary. This is the case for 
all junctions, which we believe has its origin in the fact 
that magnetization precession and possibly spatially non- 
uniform spin- wave modes of the free layer, are excited at 
currents just below this switching boundary (see Fig. [l}. 
Above the demagnetization field (~ 1.5 T), the threshold 
current increases with the applied magnetic field, consis- 
tent with Eq. [3[ Fig. [3Jb) shows the steps in V/I at 
3, 5 and 7 T. The hysteresis decreases with decreasing 
field, and vanishes at H*. H* is nearly independent of 
thickness, and is 2.5 ± 0.3 T for most junctions 15j. The 
change in DC resistance SR/R at the threshold current as 
a function of field is shown in Fig. [He). The solid line at 
2.7 % is the in-plane MR. SR/R increases with increasing 
field and asymptotically approaches the in-plane MR. In 
the macrospin model, more abrupt changes in magneti- 
zation state occur at higher magnetic field, due to the 
increasing importance of the angular dependence of the 
spin transfer torque, parameterized by 77 in our discus- 
sion. For example, for fields greater than A-kM c r/ti, the 
model predicts switching between AP and P states at 
Jt-ap an d thus the full MR. Fig. [3Jb) also shows that 
there is slight increase in the resistance at high current 
as the field increases. 

About 20 junctions with different t and lateral dimen- 
sions were measured. The critical-current densities J c as 
a function of magnetic field for all junctions were mea- 
sured and analyzed, and 4 of the junctions, representing 
different free-layer thicknesses, are plotted in Fig. 0|a). 
From the figure, it is clear that the critical current in- 
creases with thickness at fixed magnetic field. The criti- 
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FIG. 4: (a) Solid points: critical-current density vs perpen- 
dicular magnetic field of 4 of the 50 nm series junctions with 
t — 1.9, 2.8, 4.3, 5.3 nm (from left to right). Open circles: 
critical current densities of all 50 nm series junctions extrapo- 
lated to zero t. Crosses: those of 70 x 140 nm 2 junctions, (b) 
Critical current densities as a function of free layer thickness. 
Squares: 7 T. Diamonds: 3 T. Open, gray and black data 
points are the ones of 50 x 50, 50 x 100 and 70 x 140 nm 2 
junctions respectively. Solid lines: linear fits of critical cur- 
rent densities vs t of 50 nm series. Dashed lines: those of 
70 x 140 nm 2 junctions. 



cal current density at a fixed field for all junctions is plot- 
ted in Fig. [DJb) as a function of a free-layer thickness, 
with different symbols (open, gray, and black) denoting 
the three studied junction sizes, as shown in the caption. 
The critical current density increases linearly with thick- 
ness, with a slope that increases with increasing field. 
This linear dependence is similar to the results of [12| . 

One interesting result from Fig. [31(b) is that the critical 
current density J c also shows an observable dependence 
on lateral dimension. In our data, critical current densi- 
ties of 70 x 140 nm 2 junctions are generally lower than 
those of the 50 x 50 and 50 x 100 nm 2 junctions (de- 
noted the 50 nm series for ease of reference). As a result, 
50 nm series (solid lines) and 70 x 140 nm 2 junctions 
(dashed lines) were fit separately. In order to show the 
trend within the same series, 4 typical switching bound- 
aries only from the 50 nm series were plotted in Fig. Hfa). 
From Fig. gjb), the J c shift of the 70 x 140 nm 2 junc- 
tions from that of the 50 nm series is independent of the 
free-layer thickness, as the linear fits of the two different 
lateral series are parallel to each other. 

We find that the intercepts of J c vs t at different fields 
are not zero, as seen from the linear fits in Fig. Efb). 
The intercepts, denoted by J c o, of the 50 nm series are 
plotted as circles in Fig. [3Ja), which is the "switching 
boundary in the limit of zero free-layer thickness? The 
corresponding boundary of the 70 x 140 nm 2 junctions is 
plotted as the crosses in the same figure. J c q is field de- 
pendent, which shows that dJ c o/dH is nonzero. Because 



of the curvature of J c vs H at high field for thicker free 
layers, dJ c /dH becomes field dependent, and dJ c /dH vs 
t turns out to be nonlinear. Therefore, the extrapolation 
of dJc/dH to zero thickness becomes ambiguous. Inter- 
estingly, J c at a fixed field shows a linear dependence 
on thickness despite the curvature (Fig. Q|. This linear 
dependence unambiguously gives a finite intercept of J c 
at zero free-layer thickness. Note that although the cur- 
vature of the measured switching thresholds increases at 
high fields for thicker free layers, the slope of J c q vs H is 
constant from 2 to 7 T. 

J C Q is affected by the lateral dimension, since there is a 
shift between the boundaries of the two series. However, 
these two boundaries are parallel to each other, which 
shows that d J c o /dH is independent of the lateral dimen- 
sion. dJ c o/dH = (1.5 ± 0.3) x 10 7 A/cm 2 T for both 
zero free-layer thickness boundaries, while Eq. 3 predicts 
dJ c0 /dH = - n)M s t\ t=a = 0. 

The lateral dimension of the pillars affects the critical 
currents, in that a constant shift of the critical-current 
densities between the 50 nm series and 70 x 140 nm 2 was 
observed. Micromagnetic analysis of thin disks [TH, Il7j 
shows that the principle mode of the free layer is in- 
fluenced by the lateral dimension and aspect ratio. As 
the result, the mode frequency depends on the lateral 
dimension. Since the critical current is proportional to 
the mode frequency, the change in the normal mode fre- 
quency is expected to shift the intercept of J c vs H, while 
not changing dJ c o/dH. Experimentally dJ c a/dH is not 
found to depend on lateral dimension. Furthermore, in 
contrast to the prediction of Eq. [3j the threshold bound- 
aries of most junctions do not extrapolate to — 47rM c ff at 
zero current. Similarly, this may also be associated with 
the shift of the normal mode precessional frequency. 

We note that Brataas et al. has recently computed the 
critical current for the excitation of non-uniform spin- 
wave modes [l8[ . The critical current density, considering 
the effect of the non-uniform modes, was found to depend 
linearly on magnetic-layer thickness -provided the exci- 
tation of such modes is opposed by bulk Gilbert damping. 

The magnetic anisotropic field associated with the easy 
plane anisotropy is given as 47rM e ff. Based on an FMR 
study on extended Co films 47rM e ff is a function 

of Co thickness. In the thickness range studied here, 
4ttM c s changes as 1/i, and it increases only by 0.2 T 
when Co thickness decreases from 5.3 to 1.9 nm. So 
anisotropy is also a minor effect, which cannot be the 
main contribution to the non-zero dJ c o/dH. 

There are two possible explanations for the observa- 
tion of nonzero dJ c o/dH. The first is associated with 
an interface contribution to the magnetization damping. 
Recently, Tserkovnyak et al. [ll[ employed a scatter- 
ing theory approach to characterize this contribution to 
the damping for a thin ferromagnetic layer in contact 
with normal metals. They consider the spin current into 
adjacent normal layers (N) when there is magnetization 
dynamics. When such a spin current is dissipated by 
spin-flip scattering in the N layers, it generates additional 
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damping. We can estimate this additional contribution 
to the damping from the Tserkovnyak et at theory, as- 
suming that the Co layer is surrounded by perfect spin 
sinks, one of which is the fixed magnetic layer, and the 
other of which is a Pt layer, separated by 10 nm of Cu 
from the thin Co layer. The additional spin-pumping 
contribution to the damping is a' ~ 1.7 x 10~ 2 nm/t, 
and is consistent with our recent FMR studies on sim- 
ilar structures [l9|. As a result the net damping is 
a = a a + a' , where a a is the bulk damping. With an 
interface contribution to the damping, the threshold cur- 
rent goes to a finite value in the zero thickness limit. 
Quantitatively, taking rj = 0.3 and P D = 0.3, we find 
dJ c o/dH = 1.4x 10 7 A/cm 2 T, which is to be compared to 
our experimental result of 1.5 x 10 A/cm 2 T. This quan- 
titative agreement suggests interfacial damping plays a 
significant role in determining the critical current den- 
sity in spin-transfer devices. 

An alternative explanation for our results is due to 
Zhang and Levy 0, [l(| • In their model, the transverse 
component of spin decays on a length scale of Aj in the 
ferromagnet, which they find is about 3 nm for Co. As a 
result, for layers of order and less than Aj, the efficiency 
of the spin-transfer torque decreases -as the transverse 
component of angular momentum is not fully transferred 



to the thin magnetic layer. cLJt-ap / 'dH decreases to a 
t = limit of ^Aj(l-77)A/ S = 1.0 x 10 7 A/cm 2 T, with 
a = 0.01 This is also close to what we observe in 
experiment. 

In summary, the phase diagrams for spin-transfer- 
induced magnetic excitations are experimentally deter- 
mined in perpendicular magnetic field as a function of 
free-layer thickness. Based on the macrospin model (Fig. 
[1]) and the observed hysteresis, we estimate the ratio of 
the torque in the AP state to that in the P state, due to 
its angular dependence, to be approximately 2 for a con- 
stant current. Further, we have found that the critical- 
current density for spin-transfer excitations is a linear 
function of free-layer thickness that extrapolates to a fi- 
nite critical current at zero free-layer thickness. We have 
highlighted the role of the spin-pumping contribution to 
the damping, which can quantitatively explain our re- 
sults. An implication of these results is that reducing 
the thickness of the magnetic layers permits a reduction 
of the critical-current density only to a lower limit set by 
interface effects. 

We thank P. M. Levy for many useful discussions of 
this work. This research is supported by NSF-DMR- 
0405620 and by ONR N0014-02-1-0995. 
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